Renormalization group transformations on quantum states 
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We construct a general renormalization group transformation on quantum states, independent of 
any Hamiltonian dynamics of the system. We illustrate this procedure for translational invariant 
matrix product states in one dimension and show that product, GHZ, W and domain wall states are 
special cases of an emerging classification of the fixed points of this coarse-graining transformation. 
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The Renormalization Group (RG) provides a proce- 
dure to obtain an effective long distance description of 
a physical system. Following Wilson's seminal ideas 0, 
RG transformations are usually constructed in the space 
of Hamiltonians and are made of two distinct steps: first, 
a coarse-graining transformation is implemented to in- 
tegrate out short-distance information and, second, a 
rescaling of length scales and operators restores the origi- 
nal picture. This transformation is exact in the sense that 
long-distance observables remain unaltered, since they 
can be computed either with the original operators and 
Hamiltonian or with their renormalized counterparts to 
yield the same result. The exact RG transformation can 
be conveniently truncated so as to have a very powerful 
technique to retain only relevant long-distance degrees of 
freedom. 

The success of RG is ubiquous. Wilson's original idea 
has been modified such as to yield the Density Matrix 
Renormalization Group (DMRG) ^] algorithm which op- 
timizes the RG truncation; that is, the choice of rele- 
vant degrees of freedom to be retained. DMRG has been 
highly successful in describing the ground state proper- 
ties in one dimensional non-critical systems, and it can 
be understood as a variational method within the set of 
so-called matrix product states discussed below. 

^From a Quantum Information (QI) perspective, re- 
newed attention has been placed on quantum states. 
Many unexpected quantum state properties have been 
discovered and analyzed irrespectively of the dynamics 
that may produce them. It is then natural to review our 
understanding of RG focusing only on quantum states. 
This is actually the very origin of KadanofF's classical 
block spin transformation which has not been pur- 
sued on quantum states so far. The reason for the lack of 
a quantum coarse-graining analysis is related to the dif- 
ficulty of parametrizing the Hilbert space of many-body 
systems as opposed to simple Hamiltonians and to the 
complexity of dealing with wave functions in quantum 
field theory. 

Let us now introduce the general idea of the RG trans- 
formation on pure quantum states for m-party systems, 
where every local degree of freedom corresponds to a 
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Figure 1: RG evolved states can either be computed by apply- 
ing Wilsonian RG to Hamiltonians and then computing their 
ground states or by constructing quantum state RG transfor- 
mations. 



d-dimensional Hilbert space. In analogy with standard 
RG, we proceed in steps in which we: (i) merge groups 
of neighboring particles into new ones, and rescale the 
variables correspondingly; (ii) identify states which are 
equivalent under local unitary operations. This identifi- 
cation is motivated by the fact that physics at long scales 
does not depend on the choice of local basis and, as it 
will become clear below, gives rise to coarse-graining and 
irreversibility. Technically, (ii) is realized by introducing 
an equivalence relation in Hilbert space, namely = ^ if 
3Ui, . . . , Um such that lip) = J7i • • • (8) where Ut 

are local unitary transformations; that is, two states are 
equivalent if they differ by a change in the local basis. 
Thus, the RG transformation in each step can be viewed 
as a map between the resulting equivalence classes {ipo} 
(Fig. 1). In practice, we perform the RG transformation 
on a representative which is conveniently selected after 
each step. 

We describe now the above procedure in more detail 
for a ID system with translational symmetry. Given 
the representative of a class ^po, we determine the rep- 
resentative of the class in the next step, ipQ, as follows. 
We pairwise group the sites in the system and define 
a coarse-graining transformation for every pair of lo- 
cal basis states, e.g. for the sites 2j and 2j + 1, as 
\p)2j\q)2j+i — \p<l)j- This transformation yields ipo ip. 
Then we have V'o = U ^ . . . <Si U\i{j), where the (P x (P 
unitary matrix U performs the change of representative 
in the coarse-grained space. Let us emphasize that the 
freedom to take a unitary transformation in the coarse- 
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grained spaced goes beyond the onset freedom made by 
the product of unitary matrices in the original Hilbert 
spaces. The matrix U can be non-local as seen by the 2j 
and 2j + 1 sites. Some local information is now washed 
out, while preserving all the quantum correlations relat- 
ing the coarse-grained block to other ones. 

Operators also get coarse-grained along the above 
transformation. Take for instance an operator acting on 
one local Hilbert space, e.g. 02j. Expectation values 
must remain unchanged, 

{yjo\02M = mo'.m (1) 

which leads to 

o; = u{02j m2j+i)u^ (2) 

where I is the identity matrix. To complete a RG trans- 
formation we simply need to rescale distances, i.e., to 
double the lattice spacing. 

Exact RG transformations are often truncated in order 
to become of practical use. Let us, for example, consider 
the state after n coarse-graining steps, where each local 
site corresponds to 2" original spins. We are interested 
in describing long-range effects, hence in the degrees of 
freedom of the cP dimensional Hilbert space that couple 
to the outer sites (all the other information is local). In 
the case of ground states of 1-D noncritical spin chains 
for example, we know that the entropy of a block of spins 
saturates to a finite value, indicating that indeed very few 
degrees of freedom couple to the outer sites (Qj. Mathe- 
matically, this means that there exists a local isometry U 
that transforms the (P dimensional space into a much 
smaller one, without practically affecting the correlations 
between the different blocks. The whole issue of trunca- 
tion in RG then consists of coming up with an optimal 
algorithm for keeping the relevant degrees of freedom. 
Obviously, these relevant degrees of freedom will exactly 
be the ones that correspond to the largest weights in the 
reduced density operator of the block of 2" spins. This is 
very much related to the concept of the density operator 
in DMRG, although in that case only half- infinite blocks 
are considered. 

In the following we explicitly carry out the RG trans- 
formations in terms of matrix product states (MPS). In 
this representation, the whole procedure can be naturally 
implemented, and the nature of the fixed points becomes 
more transparent. Any ID translationally invariant state 
can be written in its MPS form as 

d 

IV) = J2 T^HAP^AP^...AP-)\p,,p2,...,p.m) 

(3) 

where the D x D matrices {A^} parameterize the state. 
The value of the dimension D < d"^^^ depends on the 
particular state 0. 
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Figure 2: Coarse-graining matrix product states 

The quantum coarse-graining procedure where we map 
two neighboring spins to one new block spin can be fully 
characterized in terms of the matrices A for the corre- 
sponding representatives. These matrices can be conve- 
niently chosen starting from the coarse-grained matrices 
j{{pq) ;= A^A'^ through the singular value decomposition 

min(ci^,_D^) 

(A'^''^,^ ^ {U^t'h^iV^U. (4) 
1=1 

^From this decomposition we identify the isometry U 
which selects the representative and the coarse-grained 
tensor A' 

AP ^A^ ^ X^VK (5) 

The advantage of this representative is that the Hilbert 
space corresponding to the block spins remains bounded 
above by at any step, as it is clear from the de- 
composition Q). That is, by coarse-graining we do not 
have to increase the dimension of the spins once we 
reach D^, and therefore, it is possible to perform an ex- 
act coarse-graining on finite dimensional matrix product 
states without any need of truncation! Obviously, the 
interesting question is now to classify all possible fixed 
points of this exact renormalization flow. 

Before introducing a formalism to characterize those 
fixed points, let us present some simple examples. The 
first one is provided by product states (for which D = 1, 
= 1 and = for fc = 2, . . . , d and A'^ = A*). These 
are precisely the state obtained for massive theories in 
their infrared fixed points. A significant further example 
corresponds to GHZ-like states jSi] . If we take d = 2 and 
A° = |0)(0|, A^ = |1)(1| the RG transformation reads 

A' ^ A'' = A' with J7 = |0)(00| + |1)(11|. Those states 
indeed appear as RG infrared points in e.g. the quantum 
Ising chain for vanishing external magnetic field [9| . 

Let us now construct a general formalism by which the 
fixed points can be characterized. In order to circumvent 
the arbitrariness in the choice of the local bases, we in- 
troduce the auxiliary x transfer matrix 

D 

p=i 
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where the bar indicates complex conjugation. Note that 
if we choose Aq = J^p U^A^ with U unitary, we have 



^ = ^i«® A'^= ^ (g) = 



(7) 



Conversely, the matrix E uniquely defines the matri- 
ces AP up to such a local unitary operator, and thus 
it parameterizes the equivalence class of D-dimensional 
MPS where all elements of the class are related by local 
unitary operations 11]. Furthermore, a RG step corre- 
sponds to the simple transformation E E' = E^. 
Therefore, in order to study the alluded fixed points we 
just have to characterize the class of possible operators 
{Eoo] = {liin„^oo -E"} for some E in the form Since 
we can always choose the largest eigenvalue of E to be 
equal to 1 [I4I, we just have to characterize all matrices E 
of the form ((HJ that have only eigenvalues of magnitude 1 
or (eigenvalues smaller than 1 will decay exponentially 
along repeated coarse-graining steps). 

In the generic case, the largest eigenvalue of E will 
be non-degenerate and both its left and right eigenvec- 
tors will have maximal Schmidt rank; that is, E°° = 
\^r){^l\, where the reduced density operators of $L,fl 
are rank D matrices. With the similarity transforma- 
tion fl3| we can always choose \^r) — J2iLi N) ^^'^ 
= X]fc=i -^il"); "^ith Xi > 0. From these eigenvectors 
we can directly read off the matrices A^p'^^ = yOy^\p){q\ 
with p, q = 1, . . . , D. Thus, at the fixed point, the scale- 
invariant state (representative of the equivalence class) 
can be written in terms of two auxiliary spins at each 
point, which are in an entangled state with its neigh- 
bors as shown in Fig. 2(c). This conclusion is very ap- 
pealing, and has the following consequences: 

• All connected correlation functions of the form 
(O^O^) — (O*) (O^) are exactly zero when j > i + 
this is exactly what one expects from a generic state 
by looking at it in a coarse-grained way: the corre- 
lations decay exponentially along the renormaliza- 
tion flow and become zero at the fixed point. 

• The entropy of a block of spins of length L is in- 
dependent of L and is exactly twice the entropy of 
entanglement of |$_r). 

• In terms of the picture introduced in T], coarse- 
graining turns the virtual underlying spins into real 
ones but changes the maximally entangled states 
into 

The qualitative features of the renormalization flow 
can also be easily understood. For example, an observ- 
able that only acts nontrivially on the spins in the center 
of a block converges exponentially fast to the identity 



operator times its expectation value; this was indeed ex- 
pected as observables defined in the middle of a block of 
length L with L larger than the correlation length should 
not be able to act nontrivially on the spins far away. 
On the other hand, the entropy of a block increases by 
coarse-graining (if we do not rescale the lengths). By the 
strong subadditivity of the entropy (cf.0|) this is even 
true in general for all translationally invariant states. 

As an illustration of the generic case analyzed above, 
let us consider the ground state of the AKLT Hamilto- 
nian for which A^ — {p = 1..3) are the Pauli 
matrices. Some simple algebra shows that E has eigen- 
values 3, 1, —1, —1 and hence no degeneracy in the largest 
eigenvalue. The fixed point is given by the eigenvector 
corresponding to eigenvalue 3, namely Ei ~ |/)(/| with 
/) = X]i=i I**)- Hence, it is a dimer of maximally entan- 
gled states, i.e., in the language of QI a perfect resource 
for a quantum repeater. Intriguingly, this is also the ex- 
act fixed point of the 1-D cluster state represented 
by A° = 10)((0| + (11) , A^ = |1)((0| - (1|). However, 
in this case the fixed point is already obtained after one 
coarse-graining step as the corresponding E contains a 
nilpotent Jordan block: 
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11 1 1 \ 
0000 
0000 
VI -1 -1 1/ 



E',, = J2\^^){JJ\. (8) 



The formalism developed above can also be applied 
in the case D — > 00. To illustrate that, let us con- 
sider an Ising chain with transverse magnetic field or 
an anisotropic Heisenberg antiferromagnet. In the non- 
critical regime, their correlation length ^ is finite, which 
implies that the associated operator E will have a non- 
degenerate largest eigenvalue {!/£_ oc ln(Ai/A2)). Ac- 
cordingly, the fixed point of the coarse-graining map will 
therefore consist of a dimer state. The eigenvalues of 
the reduced density operator of a half-infinite chain are 
precisely the Schmidt coefficients of the bipartite states 
making up the dimer. We can in fact determine these 
coefficients by using the results of Peschel and col. [l^ . 
who have been able to determine the mentioned reduced 
density operator exactly. For the Ising model we have 

_ re-i:r^o(2.^+iK A<i 

\»i,"2.-n^) - S g-ei:~o2j», x>i. 

Here e = 7ri^[(l - fj.'^)'^/'^]/ K{fi), with K the complete 
elliptic integral of the 1st kind, /i = min(A, 1/A), and 
rij £ {0, 1}. In the case of the noncritical Heisenberg 

model Hxxz = E^ (^'f^^'f+i + «+i + with 
A > 1, the eigenvalues are as in the Ising case with A > 1, 
but with e = arccosh(A). 

We consider now the full classification of fixed points 
for the simplest non-trivial case D = 2. By considering 
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the right eigenvectors of E corresponding to the maxi- 
mal eigenvalue, we have two possibihties [ij]: (a) one of 
them, say \^r), is entangled; (b) there exists only one 
eigenvector which corresponds to a product state. 

In the first case (a), by using the similarity transfor- 
mation we can always take \^b.) ~ X]i=i N)- Using 
the isomorphism and the full classification of trace 
preserving completely positive maps acting on a qubit 
|l4j one can prove that Eac has rank 1,4 or 2. For the 
first cases, we recover the generic case studied above or 
obtain a product state, respectively. When the rank is 2 
we obtain A° = |0)(00|; = i-e., the GHZ state 

studied above. 

The situation is more complicated in the second case 
(b), as this implies that E is not diagonalizable but has a 
Jordan-block decomposition. Some tedious but straight- 
forward algebra leads to two different possibilities for the 
fixed points. In the first case the effective Hilbert space 
is 2-dimensional and 
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exp(-i6i) 
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(10) 



Note that A^A^ 



e '^A^ and that A^ is nilpotent; 



this immediately implies that the state, when writ- 
ten in the computational basis, consists of terms like 
00 • • ■ 010 • • • 0) where at most one |1) appears. When 
6* = we recover the well-known W^-state which is 
indeed scale invariant. In the second case of Jordan fixed 
points, the spins have effective support on a Hilbert space 
of dimension 3 (but 2 when cos (a) = 1), and a possible 
decomposition is given by 



A'= , ^ ... A- = 





cos(q:) sin(/3) e 

e^*^ 
cos(a) cos(/3) 




sin(a) 



(11) 



As A^A^ (X A^^A^A^ oc A^.A'^A^ oc A^.A^A^ = 
0,^2^1 = 0,^^^" = 0, the state will be a superposi- 
tion of terms of the form |00 • ■ • 0122 • • ■ 2). Therefore 
these scale-invariant states represent linear combination 
of domain walls. 

For the case D — 2 this completes the classification of 
all fixed points, which remarkably correspond to almost 
all well-studied multipartite states encountered in Ql. We 
note that in the thermodynamic limit case (b) in the clas- 
sification is redundant in the sense that W and product 
state as well as domain wall and GHZ state become then 
locally indistinguishable. A distinction can, however, be 
relevant from a QI perspective, where it is more natural 
to carry out only a finite number of coarse-graining steps 
which enlarge the region of local accessibility only to a 
physically reasonable vicinity. In this context the system 
to be transformed can thus be finite. 

For D > 2 a coarse classification of the fixed points can 



be given by their possible decompositions into ergodic 
states and their periodic components discussed in 

Finally, let us mention that we have concentrated here 
on ID systems, since in that case MPS give a very useful 
description. Fortunately, the analogue of MPS in higher 
dimensions has been recently put forward |l8j | , which may 
allow us to extend the analysis introduced here to two 
and three spatial dimensions. 

Let us conclude with a comment on irreversibility of 
RG flows. Our RG transformation defined on finite ma- 
trix product states incorporates unitarity in a natural 
way at variance with other approaches. Moreover, the 
eigenvalues of E^ are smaller than the ones of E. RG 
flows will inevitably eliminate all the eigenvalues smaller 
than 1, which can be phrased as irrelevant pieces of the 
state. A careful treatment of irreversibility of RG flows 
should incorporate the infinite dimensional case, associ- 
ated to conformal field theories, and the discussion of 
degeneracy in the final relevant Hilbert space. 
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